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Abstract 

Inclusive charged hadron production, e"'"e~ e"'"e~ h^ X, is studied using 414 
pb~^ of data collected at LEP with the L3 detector at centre-of-mass energies be- 
tween 189 and 202 GeV. Single particle inclusive differential cross sections are mea- 
sured as a function of the particle transverse momentum, pt, and pseudo-rapidity, 
7]. For Pt < 1.5 GeV, the data are well described by an exponential, typical of 
soft hadronic processes. For higher pt, the onset of perturbative QCD processes is 
observed. The tt^ production cross section for pt > 5 GeV is much higher than the 
NLO QCD predictions. 
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1 Introduction 



Two-photon collisions are the main source of hadron production in the high-energy regime of 
LEP via the process e^e~ — >■ e^e~7*7* e^e~hadrons. In the Vector Dominance Model 
(VDM), each photon can transform into a vector meson with the same quantum numbers, thus 
initiating a strong interaction process with characteristics similar to hadron-hadron interactions. 
This process dominates in the "soft" interaction region, where hadrons are produced with a 
low transverse momentum, pt, with respect to the beam direction. Hadrons with high pt are 
produced by the direct QED process 7*7* —>■ qq or by QCD processes originating from the 
partonic content of the photon. QCD calculations are available for single particle inclusive 
production in two-photon interactions at next-to- leading order (NLO) precision [HEl- 

The L3 collaboration recently published results on inclusive vr*^ and Kg production [3] . The 
7r° differential cross section measured as a function of pt exhibits a clear excess over QCD 
calculations. A comparison of these results with other single particle inclusive production 
at high Pt is therefore important. In this Letter, the inclusive charged hadron production is 
studied for a centre-of-mass energy of the two interacting photons, W^-y, greater than 5 GeV. 
The hadrons are measured in the transverse momentum range 0.4 GeV < < 20 GeV and in 
the pseudo-rapidity^^ interval |?7| < 1. The contributions from vr^ and are also derived. 

The data used for this analysis were collected by the L3 detector |3] at centre-of-mass 
energies y/s = 189 — 202 GeV, with a luminosity weighted average value of y/s = 194 GeV, 
for an integrated luminosity of 414 pb~^. Results on inclusive charged hadron production for a 
smaller data sample at lower ^/s were previously reported jS]. 

The process e"'"e~ e^e~ hadrons is modelled with the PYTHIA [H] event generator for an 
event sample three times larger than the data. In this generator, each photon can interact as a 
point-like particle, as a vector meson or as a resolved photon, leading to six classes of events. 
The fragmentation is simulated with JETSET. Predictions from the PHOJET Monte Carlo 
program [2j are also compared with the data. The following Monte Carlo generators are used 
to simulate the relevant background processes: KK2f jH] for e+e"— *• qq(7); KORALZ [HI for 
e+e"^ r+r-(7); KORALW [10] for e+e"^ W+W" and DIAG36 [H] for e+e"^ e+e-r+r". 
Events are simulated in the L3 detector using the GEANT and GHEISHA ^3] programs 
and passed through the same reconstruction program as the data. Time dependent detector 
inefficiencies, as monitored during each data taking period, are also simulated. 

2 Event and charged hadron selection 

Two-photon events are collected predominantly by the track triggers ^1] with a low pt threshold 
of about 150 MeV. The selection of e"'"e~ e^e~ hadrons events fSj consists of: 

• A multiplicity cut. To select hadronic final states, at least six objects must be detected, 
where an object can be a track or a calorimetric cluster with no associated track. 

• Energy cuts. The total energy deposited in the calorimeters must be less than 0.4 a/s, 
in order to exclude e"'"e~ annihilation events. The total energy in the electromagnetic 
calorimeter is required to be greater than 500 MeV, to suppress beam-gas and beam-wall 
backgrounds. 

^^?7 = — In tan(6'/2), where 9 is the polar angle of the particle relative to the beam axis. 
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• An anti-tag condition. Events with a cluster in the luminosity monitor with an energy 
greater than 30 GeV and an electromagnetic shower shape are excluded. 

• A mass cut. The visible mass of the event must be greater than 5 GeV. 

About 2 million hadronic events are selected by these criteria. The overall background level 
is less than 1% and is mainly due to the e"'"e~ — > qq (7) and e+e~— >• e'^e~T~^r^ processes. 

Charged hadrons are measured with high quality tracks in the inner tracking detector. 
These tracks have a transverse momentum greater than 400 MeV and a distance of closest 
approach to the primary vertex in the transverse plane less than 4 mm. The number of hits 
must be greater than 80% of that expected from the track length. Tracks are analysed in 
the |?7| < 1 and pt < 20 GeV range where the detector resolution is optimal. A resolution 
apjpt ~ 0.015( GeV"^) x pt is achieved. 

3 Differential cross section 

The differential cross sections of inclusive charged hadron production as a function of pt are 
measured for an effective mass of the 77 system W^^ > 5 GeV, with a mean value of {W-y^) ~ 
30 GeV, a photon virtuality <S GeV^ and an average photon virtuality (Q^) ~ 0.2 GeV^. 
This phase space is defined by cuts at the Monte Carlo generator level. Results are presented 
in 12 Pt bins between 0.4 and 20 GeV. 

The distribution of the detected charged hadrons in these pt bins is presented in Figure C^. 
The background remains very low over the whole pt range. Events from the e"'"e~— > e'^e~T~^T~ 
process dominate the background at low pt while annihilation events dominate it at high pt. To 
measure the cross section, the background is subtracted bin-by-bin and the data are corrected 
for the selection efficiency, including acceptance, calculated bin-by-bin with PYTHIA. This 
selection efficiency varies from 62% to 84%. At low pt, the efficiency decreases due to the effect 
of the mass and energy cuts. At high pt, it decreases because of the multiplicity cut, since high 
Pt particles are mainly produced in low multiplicity events. 

The level 1 trigger efficiency is obtained by comparing the number of events accepted by 
the independent track and calorimetric energy jTHj triggers. It varies from 95% to 98%. The 
efficiency of higher level triggers is about 95% and is measured using prescaled events. The 
overall efficiency, taking into account selection and trigger efficiencies is given in Table ^ 

Sources of systematic uncertainties on the cross section measurements are the trigger effi- 
ciency estimation, the background subtraction, the selection procedure and the Monte Carlo 
modeling. Their contributions are shown in Table El The uncertainty on the trigger efficiency 
and on the background subtraction are of a statistical nature. The uncertainty due to the 
selection procedure is evaluated by repeating the analysis with different selection criteria: the 
multiplicity cut is moved from 5 to 7 objects, the energy cut is moved to 0.35 y/s and the 
number of hits of the tracks is moved to 70% of that expected. The sum in quadrature of the 
differences between these and the reference results is listed in Table El Varying other criteria 
give negligible contributions. To evaluate the uncertainty on the Monte Carlo modeling, the 
selection efficiency is determined using only one of the PYTHIA subprocesses: VDM-VDM, 
direct-direct or resolved-resolved. The systematic uncertainty is assigned as the average differ- 
ence between these values and the reference Monte Carlo. The larger contribution comes from 
the difference between direct and other processes. 

The differential cross section of charged hadron production as a function of pt is presented 
in Figure and in Table ^ The migration due to the pt resolution does not affect these 
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results. This was verified by performing a one-step Bayesian unfolding ^\ of the track pt 
distribution which give results compatible, within errors, with those obtained using the bin-by- 
bin correction. 

The steep decrease of da/dpt in the range 0.4 < pt < 1.5 GeV is described by an expo- 
nential of the form Aexp{—pt/{pt)) with a mean value of (pt) — 232 MeV. This behaviour 
is characteristic of hadrons produced by soft interactions and is similar to that obtained in 
hadron-hadron and photon-hadron collisions P^^. At higher pt the differential cross section is 
better represented by power law functions Ap^^ , as expected by the onset of QCD processes. 
For 1.5 GeV < Pt < 5 GeV, B ~ 4.2 and for 5 GeV < pt < 20 GeV, B ~ 2.6. The results of 
the fits are drawn on Figure^ where the data are also compared to Monte Carlo predictions. 
PYTHIA is slightly above the data, whereas PHOJET is too low by more than one order of 
magnitude. These results are consistent with our findings in inclusive 7r° production jS] . 

4 Charged pions and charged kaons 

Assuming the fragmentation function implemented in JETSET are correct, the vr^ and the 
inclusive cross sections are extracted from the charged hadron cross section. Their ratios 
relative to charged hadrons are estimated bin-by-bin from Monte Carlo. Above 5 GeV, they 
are almost constant. Their uncertainty is calculated in the same way as the uncertainty on the 
Monte Carlo modeling of the selection efficiency, by using different subprocesses in PYTHIA. 
This gives an additional systematic uncertainty of from 2% to 12% for pions and from 14% to 
24% for kaons. 

The differential cross sections for vr^ and production as a function of pt are presented 
in Figure El and in Table El The vr^ data are compared to the previous 7r° data jS] scaled up by 
a factor 4: a factor 2 to correct for the \ri\ < 0.5 interval used for the vr" measurement and a 
factor 2 to take into account the isospin symmetry. A good agreement is found between these 
two measurements as shown in Figure El The data are compared to the previous results of 
Kg data scaled up by a factor 4/3: a factor 2/3 to correct for the |?7| < 1.5 interval of the 
Kg measurement and a factor 2 to take into account unobserved Kl decays. Good agreement is 
found between these two measurements as shown in Figure El These agreements show a good 
consistency with the data of the fragmentation functions as implemented in JETSET. 

The differential cross section of tt^ production as a function of |?7| for pt > 1 GeV is shown 
in Figure El and in Table O The cross section is almost constant in this r] range. It agrees well 
with the 7r° measurement 0. For different pt cuts, Monte Carlo and QCD predictions describe 
well the uniform rj distribution, while the agreement in the absolute rate depends on the pt 
range considered. 

In Figure Et. the data are compared to analytical NLO QCD predictions p^l2]. For this 
calculation, the flux of quasi-real photons is obtained using the Equivalent Photon Approxi- 
mation |2IIj, taking into account both transverse and longitudinal virtual photons. The inter- 
acting particles can be point-like photons or partons from the 7 ^ qq process, which evolve 
into quarks and gluons. The NLO parton density functions of Reference are used and all 
elementary 2 — 2 and 2 — 3 processes are considered. New NLO fragmentation functions j22| 
are used. The renormalization, factorisation and fragmentation scales are taken to be equal: 
H = M = Mp = ^Pt |2|, with ^ = 1 for the central value. The scale uncertainty in the NLO 
calculation is estimated by varying the value of ^ from 0.5 to 2.0. The agreement with the data 
is poor in the high-pf range for any choice of scale. 

To test NLO QCD calculations in regions where non-perturbative subprocesses are better 



4 



suppressed, we have also measured differential cross sections of vr production for W^^ > 10, 
30 and 50 GeV. The results are shown in Table |3 and Figure Eb- The discrepancy between the 
calculations and data at high pt is not significantly reduced by these stringent more W^.y cuts. 

Similar calculations were previously compared to 7p reactions at HERA up to a pt of 
12 GeV and to pp collisions up to a pt of 20 GeV. Good agreement was found |23j. In the 77 
channel, an excess of data with respect to NLO QCD was observed in tagged events at PETRA 
experiments |2]. No discrepancy is observed with the OPAL data which explore a pt range 
up to 10 GeV. In this range, our data and the OPAL ones are well in agreement within the 
quoted uncertainties. A discrepancy with NLO QCD is revealed by our data which extend the 
measurement to higher pt values. 
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(10.9 


± 


0.1 ± 


1.5) 


xlO-'^ 


0.8 - 1.0 


0.88 


67.7 ± 


6.0 


(48.0 


± 


0.1 ± 


5.9) 


xlO^ 


1.0 - 1.5 


1.14 


72.4 ± 


4.8 


(14.1 


± 


0.1 ± 


1.4) 


xlO^ 


1.5 - 2.0 


1.68 


77.4 ± 


3.7 


(28.5 


± 


0.1 ± 


2.2) 


xlO 


2.0 - 3.0 


2.31 


77.2 ± 


4.3 


(60.9 


± 


0.5 ± 


4.4) 




3.0 - 4.0 


3.36 


75.0 ± 


5.2 


(13.1 


± 


0.2 ± 


1.0) 




4.0 - 5.0 


4.39 


69.5 ± 


5.8 


(48.7 


± 


1.3 ± 


4.2) 


xlQ-^ 


5.0 - 7.5 


5.79 


68.1 ± 


6.8 


(15.3 


± 


0.4 ± 


1.6) 


xlO^i 


7.5 -10.0 


8.46 


65.2 ± 


8.7 


(50.9 


± 


2.5 ± 


7.0) 


xlO-2 


10.0 -15.0 


11.98 


61.9 ± 11.0 


(21.0 


± 


1.2 ± 


3.8) 


XlO-2 


15.0 -20.0 


17.36 


59.8 ± 14.7 


(97.1 


± 


8.4 ± 24.3) 


xlO-3 



Table 1: Transverse momentum range and average value from the data with the corresponding 
overall efficiency and differential cross section for inclusive charged hadron production for W^y > 
5 GeV and [t^I < 1. The uncertainty on the efficiency is systematic. The first uncertainty on 
the cross section is statistical and the second systematic. 



Pt 


trigger 


background 


selection 


Monte Carlo 


[GcV] 


ofiicioncy ['X] 


sul)tractiou ['X] 


procedure [*X] 


modeling [%] 


0.4 ^ 0.6 


0.1 


< 0.1 


10.1 


12.4 


0.6 - 0.8 


0.1 


< 0.1 


9.2 


10.6 


0.8 - 1.0 


0.2 


< 0.1 


8.4 


8.9 


1.0 - 1.5 


0.2 


< 0.1 


7.5 


6.6 


1.5 - 2.0 


0.4 


< 0.1 


5.9 


4.8 


2.0 - 3.0 


0.5 


0.1 


4.5 


5.6 


3.0 - 4.0 


0.9 


0.4 


3.1 


6.9 


4.0 - 5.0 


1.2 


0.9 


2.3 


8.2 


5.0 - 7.5 


1.2 


1.3 


1.6 


10.0 


7.5 -10.0 


1.2 


2.5 


1.2 


13.3 


10.0 -15.0 


1.2 


3.1 


1.1 


17.8 


15.0 -20.0 


1.2 


4.2 


1.1 


24.6 



Table 2: Systematic uncertainty on the charged hadron cross section due to trigger efficiency, 
background subtraction, selection procedure and Monte Carlo modeling. 
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{pt) 


da/dpt for pions 


da/dpt for kaons 


[GeV] 






[pb/GeV] 








[pb/GeV] 




0.48 


(20.3 


± 


0.1 ± 


3.3) 


xlO^ 


(23.7 


± 


0.1 ± 5.1) 


xlO^ 


0.68 


(88.1 


± 


0.1 ± 12.5) 


xlO^ 


(15.3 


± 


0.1 ± 3.1) 


xlO^ 


0.88 


(36.5 


± 


0.1 ± 


4.5) 


xlO^ 


(80.7 


± 


0.2 ± 15.3) 


XlO 


1.14 


(10.2 


± 


0.1 ± 


1.0) 


xlO^ 


(26.4 


± 


0.1 ± 4.7) 


xlO 


1.68 


(20.5 


± 


0.1 ± 


1.6) 


xlO 


(54.3 


± 


0.4 ± 9.2) 




2.31 


(44.7 


± 


0.3 ± 


3.5) 




(10.9 


± 


0.1 ± 1.9) 




3.36 


(10.0 


± 


0.2 ± 


0.9) 




(21.3 


± 


0.4 ± 3.8) 


xlO-i 


4.39 


(37.8 


± 


1.0 ± 


3.6) 


xlO-i 


(71.3 


± 


2.4 ± 13.4) 


xl0~2 


5.79 


(12.3 


± 


0.4 ± 


1.4) 


xlO-^ 


(20.6 


± 


0.8 ± 4.2) 


xlO-2 


8.46 


(41.0 


± 


2.1 ± 


6.2) 


xlO-2 


(62.9 


± 


4.1 ± 14.8) 


xlO-3 


11.98 


(16.9 


± 


1.0 ± 


3.4) 


xlO-2 


(27.7 


± 


2.1 ± 7.7) 


xlO--'^ 


17.36 


(81.3 


± 


7.1 ± 22.5) 


xlO-3 


(10.7 


± 


1.4 ± 3.7) 


xlO-3 



Table 3: Differential cross section as a function of pt for inclusive tt and K production for 
W^^ > 5 GeV and |7^| < 1. The first uncertainty on the cross section is statistical and the 
second systematic. 



(Pt) 


da/dpt [pb/GeV] 


da/dpt [pb/GeV] 


da/dpt [pb/GeV] 


[GeV] 


W^^ > 10 GeV 


W^^ > 30 GeV 


W^^ > 50 GeV 


0.48 


(13.7 ± 0.1 ± 


3.2) 


xlO=^ 


(56.2 ± 0.2 ± 23.8) 


xlO^ 


(30.1 ± 0.2 ± 17.0) 


xlO^ 


0.68 


(60.6 ± 0.1 ± 


12.5) 


xlO^ 


(25.3 ± 0.1 ± 9.9) 


xlO^ 


(13.8 ± 0.1 ± 7.5) 


xlO^ 


0.88 


(25.7 ± 0.1 ± 


4.5) 


xlO^ 


(10.9 ± 0.1 ± 3.9) 


xlO^ 


(59.3 ± 0.6 ± 31.2) 


xlO 


1.14 


(74.3 ± 0.2 ± 


10.3) 


XlO 


(33.5 ± 0.2 ± 10.7) 


XlO 


(18.9 ± 0.2 ± 9.5) 


xlO 


1.68 


(15.6 ± 0.1 ± 


1.0) 


xlO 


(73.5 ± 0.8 ± 8.9) 




(39.4 ± 0.7 ± 9.4) 




2.31 


(36.4 ± 0.3 ± 


2.2) 




(17.2 ± 0.3 ± 2.0) 




(96.1 ± 2.3 ± 23.4) 


xlO-^ 


3.36 


(88.8 ± 1.5 ± 


5.8) 


xlO-^ 


(40.5 ± 1.1 ± 4.6) 


xlO-^ 


(25.5 ± 1.1 ± 6.5) 


xlO-i 


4.39 


(35.7 ± 1.0 ± 


2.9) 




(18.9 ± 0.8 ± 2.2) 


xlO-i 


(96.5 ± 6.0 ± 26.0) 


XlO-2 


5.79 


(11.7 ± 0.4 ± 


1.2) 


xlO-i 


(60.4 ± 2.5 ± 7.2) 


XlO-2 


(36.5 ± 2.1 ± 10.5) 


xlO-2 


8.46 


(37.8 ± 2.0 ± 


5.6) 


xlO-2 


(22.1 ± 1.5 ± 2.9) 


xlO-2 


(13.1 ± 1.2 ± 4.4) 


xlO-2 


11.98 


(16.4 ± 1.0 ± 


3.4) 


xlO-2 


(12.7 ± 0.9 ± 1.8) 


xlO-2 


(84.1 ± 7.4 ± 31.5) 


xlO-=^ 


17.36 


(78.9 ± 7.0 ± 


23.9) 


xlO-=^ 


(60.0 ± 6.3 ± 10.1) 


xlO-^ 


(61.3 ± 7.5 ± 27.2) 


xlO-^ 



Table 4: Differential cross section as a function of pt for inclusive n"^ production for |?7| < 1 
and different W^^ cuts. The first uncertainty on the cross section is statistical and the second 
systematic. 
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V 




da/d\r)\ [pb] 


0.0 




0.2 


638 ± 3 ± 80 


0.2 




0.4 


677 ± 3 ± 84 


0.4 




0.6 


693 ± 4 ± 86 


0.6 




0.8 


719 ± 4 ± 90 


0.8 




1.0 


687 ± 4 ± 86 



Table 5: Differential cross section as a function of |?7| for inclusive tt^ production for W^^ > 5 
GeV and Pt > 1 GeV. The first uncertainty on the cross section is statistical and the second 
systematic. 
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Pt [GeV] [GeV] 

Figure 1: a) Number of selected tracks per GeV in each pt bin and main sources of background, b) Inclusive charged hadron differential 
cross section da/dpt fitted with an exponential and power-law functions. Monte Carlo predictions are also presented. Statistical and 
systematic uncertainties are shown. The average pt value of each bin, (p^) , is used. 



■D 



10 



^ 10 

CD 

o 

^ 10 

Q. 
■D 



10 



10 



L3 



1 



« 



+ 



• Data n 
o Data 71^ 
■ Data K 
□ Data K° 



y 6 



a) 



Pt [GeV] 1 



CO 



10 



8 








Pt [GeV] 1 



Figure 2: a) Differential cross section da/dpt for inclusive pion and kaon production. The tt^ data are compared to the inclusive 7r° 
measurement [3] scaled by a factor 4. The K='= data are compared to the inclusive Kg measurement [3] scaled by a factor 4/3. b) Cross 
section ratios for pions and kaons. Good agreement is found with the expected values (horizontal lines). Statistical and systematic 
uncertainties are shown. The average pt value of each bin, (p^) , is used. 
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Figure 3: Inclusive 7r=^ differential cross section da/d\r]\ for pt> 1 GeV compared to the inclusive 
7r° measurement [3] scaled by a factor 2 and two Monte Carlo predictions. Statistical and 
systematic uncertainties are shown. 
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Figure 4: a) Inclusive n differential cross section da/dpt compared to NLO QCD calculations [19] for WC^^ > 5 GeV. The dashed- 
dotted line corresponds to the direct subprocess. The dashed lines represent the scale uncertainty of the calculations, b) Inclusive tt^ 
differential cross section da/dpt with different W^^ cuts. The average pt value of each bin, {pt), is used. 



